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! Abstract 

Recent results from Super Kamiokande suggest j/„ — v T mixing and hence lepton 
0^ . flavor violation. In supersymmetric models, this flavor violation may have implications 

for the pattern of slepton masses and mixings. Possible signals for this mixing in the 
1 decays of sleptons produced at the LHC are discussed. The sensitivity expected is 

. compared to that of rare decays such as r — ► /xy. 
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1 Introduction 



Recent results from Super Kamiokande and earlier measurements from Soudan-2 and 
1MB P| on the relative fluxes of electron and muon type neutrinos produced by the interaction 
of cosmic rays in the atmosphere show a deficit of muon-type neutrinos. These results are 
consistent with the oscillation of into v T provided that the two mass eigenstates are fully 
mixed, sin 2 26 m 1, and their mass difference is of order Am 2 ~ 10 _3 eV 2 . Models that 
can accommodate this mixing have lepton flavor violation (LFV) and other lepton number 
^ ■ violating processes which are likely to be observable, such as r — > /iy and r — > It 
has been suggested j|, |5j, U that supersymmetric models can naturally accommodate the 
flavor violation. 

If there is mixing of neutrinos, then one would in general expect mixing of sleptons and 
sneutrinos in SUSY. In general, SUSY models can have significant flavor mixing. The slepton 
mass matrix is a 6 x 6 matrix involving the three flavors of left and right sleptons and has 
the form 

p* ( M 2\MNp _ (p* z* \ ( M Lij M LRil \ f£ L j 

^mh t»j-{hi ^U^i Mi' kl )[~e Rl . 

where M, N = L,R label chirality and i, j,k,l = 1,2,3 are generational indices. M\ and M\ 
are the supersymmetry breaking, chirality conserving, slepton masses. The chirality breaking 
terms, M\ R - fc , contain a flavor diagonal piece mifi tan (3, where mi is the corresponding 
charged lepton mass, fi is the supersymmetry conserving Higgs mass parameter and tan/? 
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is the usual ratio of Higgs vacuum expectation values. Additional supersymmetry breaking 
terms (A terms) also appear in the chirality mixing term. These cannot be flavor diagonal 
if, as we require, the neutrino mass matrix is not diagonal in the flavor basis. However, these 
terms are expected to be proportional to charged lepton masses. We will therefore make the 
simplifying assumption that these terms are relevant only for the third generation. In order 
to accommodate the into v T mixing, there must be significant mixing in the stau, smuon 
sector ||. 

In the MSSM (minimal supersymmetric standard model) there is no lepton number viola- 
tion. It can be incorporated by adding right handed neutrinos (iVj) with couplings to the left 
handed lepton doublets Lj of the form XijNiLjH where if is a Higgs doublet. Then radia- 
tive corrections at one loop will induce lepton number violating terms in the mass matrices 
for the left sleptons. There can also be similar soft supersymmetry breaking A terms. We 
therefore assume the following form for the slepton mass matrix in the e^, /iL, t^, Sr, &r,tr 
basis IBfl : 
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where D L = -~(2M^ - Mf) cos2/3, D R = (M^ - Mf) cos 2/5, and A T = A T — fi tan/3. 
Minimal supersymmetric models such as SUGRA are recovered by neglecting M 2 T . The 
Super Kamiokande results suggest 
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It has been suggested that this large flavor violation in the slepton sector could be ob- 
served indirectly via the decay r — ► /ry ]4], [|. By contrast, this paper considers direct signals 
from the observation of lepton flavor non-conservation in the SUSY sector at the LHC We 
will demonstrate that, in certain cases, the direct observation is more sensitive than the 
flavor violating tau decay. Section |2] defines a sample point based on the minimal SUGRA 
model that has been studied in detail, and Section [5] explains how sensitivity to r — \x flavor 
violation can be observed at the LHC for this model. Section |] compares the sensitivity for 
this direct search with that from the rare decay r — > fi'y. Section [5] presents results indicating 
the range of parameters in the minimal SUGRA model for which such direct searches might 
be feasible. Finally, Section || presents a summary and conclusions. 



2 Example — Model Parameters 

We consider a minimal SUGRA 0, | point with m = 100 GeV, m 1/2 = 300 GeV, A = 
300 GeV, tan (3 = 10, and sgn /i = +. We use the ISAJET || implementation of the SUGRA 
model. This point is identical to one of those studied in detail [0, [H] (SUGRA Point 5) 
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except for the value of tan (3. The larger tan (3 gives a light Higgs mass of 113 GeV, consistent 
with the current LEP limit [|TJ]. The masses are calculated with ISA JET 7.49, which assumes 



5 = 0. A modification was made to make available the complete slepton mass matrix, giving 
access to all the soft mass parameters at the electroweak scale viz.: 

M L = 236.0 GeV, M R = 153.3 GeV, /x = 360.2 GeV 
M TL = 235.2 GeV, M TR = 150.7 GeV, A T = 94.9 GeV 

The mass eigenvalues for 5 = are 

m ih = 231.76, m iR = 156.26, m fl = 145.36, m? 2 = 232.96 GeV. 

(fi is the lighter of the two mass eigenstates resulting from the mixing of ?l and tr.) Note 
that in SUGRA type models the £r are generally lighter than the t^. The dominant final 
states for x\ decay are x\Z° , e^e and /x,R/i, and fjT which have branching ratios of 6.9%, 
12.8% (for both e and /x), and 66.1% respectively. The Xi branching ratios are 45.3% for 
XiW and 54.4% for f x v T . 

Direct Drell-Yan production of sleptons at the LHC has a small cross section (270 fb 
for all sleptons in the case considered) and is difficult to separate from SUSY backgrounds. 
In contrast the total SUSY rate, made up largely of squark and gluino production is 100 
times larger. Gluinos (of mass 706 GeV) decay to squarks (average mass 620 GeV) with 
an almost 100% branching ratio and left-squarks have a 30% branching ratio to x\ (°f mass 
212 GeV). Hence the most copious source of observed dileptons from sleptons in this case 
is x\ ~~ * where the X2 arises in squark decay. This will always be the case provided that 
the decay is allowed. If only £ R is produced in this way (in the case studied £l is heavier 
than x°)) then the main effect of 5 7^ is to add a component of /x^ in fi so that X2 ~ > ^i/ 1 
and fi — > fiXi- The x.i mass is 116 GeV. This decay gives a /i ± r =F signal and produces an 
asymmetry between /x ± r :F and e ± r T final states. 

The slepton mass matrix was rediagonalized using the soft mass parameters given above 
plus the mixing term S, inducing mixings among slepton flavors and hence flavor-violating 
decays such as X2 ~^ an d ? — ► XiA*- The resulting branching ratios for x° ~" > x\^ T an d 
for X2 ~* Xi/^A* through an intermediate f\ with at least one lepton-flavor violating decay are 
shown in Figure |I[ The latter decay is uninteresting both because it is small and because 
it cannot be separated from the X2 ~ * ~~ * TT Xi fi na l state for which both taus decay to 
muons. We therefore concentrate on the /xt decay followed by a hadronic tau decay. Note 
that the mass shifts in the lightest states are very small. For 5 = 0.1, the masses of the 
lightest states change by ~ 100 MeV; the heavier states are more strongly affected, their 
masses become 221 and 243 GeV; the mass of ei is, of course, unchanged at 231 GeV. 

A sample of 600000 SUSY events was generated with ISAJET 7.49, equivalent to an 
integrated luminosity of ~ 25 fb -1 . The detector response is parameterized by Gaussian 
resolutions characteristic of the ATLAS [0 detector without any non-Gaussian tails. In the 
central region of rapidity we take separate resolutions for the electromagnetic (EMCAL) and 
hadronic (HCAL) calorimeters, while the forward region uses a common resolution: 



EMCAL 10%/V£©1%,|77| <3 
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Figure 1: Branching ratios for ~^ Xi/ iT (solid) and for \% ~^ Xi/-*/-* through an intermediate 
f\ (dashed) as a function of the mixing parameter S. 



HCAL 50%/VE © 3%, \r]\ <3 
FCAL 100%/VE © 7%, \t]\ >3 

All energy and momenta are measured in GeV. A uniform segmentation Arj = A<f) = 0.1 is 



used with no transverse or longitudinal shower spreading. Both ATLAS |T3[ and CMS |L4 
have finer segmentation over most of the rapidity range. An oversimplified parameterization 
of the muon momentum resolution of the ATLAS detector including both the inner tracker 
and the muon system measurements is used, viz 



5p T /p T = ^0.0162 + (O.OOllpr) 2 

In the case of electrons we take a momentum resolution obtained by combining the electro- 
magnetic calorimeter resolution above with a tracking resolution of the form 




Sp T /p T = 1 + ' v /(0.0004p T ) 2 + 0.0001 



This provides a slight improvement over the calorimeter alone. Missing transverse energy is 
calculated by taking the magnitude of the vector sum of the transverse energy deposited in 
the calorimeter cells. Jets are found using GET JET |J with a simple fixed-cone algorithm. 
The jet multiplicity in SUSY events is rather large, so we have used a cone size of 



R= ^(At]) 2 + (A0) 2 = 0.4 

unless otherwise stated. Jets are required to have at least px > 20 GeV; more stringent 
cuts are often used. All leptons are required to be isolated and have some minimum p? 
and \q\ < 2.5. The isolation requirement is that no more than 10 GeV of additional be 
present in a cone of radius R = 0.2 around the lepton to reject leptons from 6-jets and c-jets. 
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Figure 2: Dilepton mass distribution: OS (opposite sign) SF (same flavor) signal (solid), 
OSOF signal (dashed), and Standard Model OSSF (shaded). 

3 Final States Involving Taus 

SUSY production at the LHC is dominated by the production of squarks and gluinos, which 
decay via cascades to the lightest SUSY particle Xi, which escapes the detector. Events were 
therefore selected by requiring multiple jets plus large $t'- 

• > 4 jets with p Tjl > 100 GeV and p T ,2,s,4 > 50 GeV; 

• M eff = $ T + Pt,i + Pt,2 + Pt,3 + Pt,a > 800 GeV; 

• $ T > 0.2M efT ; 

The same analysis was applied to the Standard Model background. ISAJET 7.44 was used 
to generate 250000 events each of it, Wj, Zj, and WW and 2500000 QCD jets, divided 
among five pt bins covering 50 < px < 2400 GeV. The cut on M e g was chosen to make this 
background fairly small but still visible. These background samples represent much smaller 
integrated luminosity that the signal sample, so the resulting statistical fluctuations are large 
when they are scaled to the appropriate integrated luminosity. 

For events that pass this selection, an additional requirement that there be two isolated 
leptons with pt > 10 GeV and \rj\ < 2.5 is made. The dilepton mass distribution for opposite 
sign, same flavor pairs is shown in Figure [| and shows the characteristic endpoint for slepton 
mediated decay. The end-point occurs at 



A small Z° peak is also visible, primarily from the 6.9% branching ratio for %° ~~* XiZ°. 
The events in the opposite sign, opposite flavor sample are due to decays of the type x° — > 



M, 



max 



(M|- Mf)(Mf - Mj, 
Mf 



95.1 GeV. 
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Figure 3: Visible r h r h mass for r^r h signal (solid), t^t^ signal (dashed), and Standard 
Model T^Th (shaded). Only hadronic tau decays are included. 
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Figure 4: Visible mass distributions for i^r^ signal (solid), P^t^ signal (dashed), [i 



,± T =F 



from LFV decays with BR = 10% (dot-dashed), and Standard Model 
hadronic tau decays are included. 



(shaded). Only 



f\T — > t + t~x1 — > e + n~v e V^v T V T and from events with leptons arising from xf decay which 
are not subject to the kinematic constraint. 

Selection of hadronic r decays was based on the actual hadronic tau decays plus the 
reconstruction efficiency and jet rejection based on the full simulation analysis done for 
an earlier study (T^|. Pile-up effects have been ignored, so the results are applicable only at 
low luminosity, ~ 10 33 cm _2 s _1 . For each hadronic r decay with > pT,mm — 20 GeV and 
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Figure 5: Visible mass distributions for E Tj[ — I (solid) and \i from LFV decays 
with BR = 10% (dot-dashed). All the Standard Model backgrounds considered here should 
cancel in this plot up to statistical fluctuations. 



7] < 2.5, the direction was found. A matching jet with rj < 2.5 was then sought with AR < 0.4 
and pt, t > 0.8;>rjet- If such a jet was found, it was tagged as a r with a probability of 66%; 
its charge was assigned correctly with a probability of 92%. Jets not so tagged were assigned 
as r's using an approximate parameterization of the jet/r rejection shown in Figure 9-31 of 



Ref. 15, 



5 

R = (o.971pf - 49 GeV) 3 



where e T is the hadronic r efficiency and px is in GeV. We choose e T = 66%. The efficiency 
and the results using the parameterization of fl6j are not very different. Note that we are 



actually using the hadronic jet resolution for all r decays here, not the improved resolution 



for a subset of decays as in the earlier analysis [17]] which was focussed on the invariant mass 
of a pair of taus, both decaying hadronically. Since this analysis relies mainly on counting 
events, the resolution on the measurement of the hadronic r decay products is not crucial. 

The visible (hadronic) rr mass for hadronic r decays (hereafter denoted by 77J is shown 
in Figure 0. Here events were required to have two r h candidates in addition to the cuts 
listed above. Note that one or both of these candidates could be the same jets that were used 
in the initial event selection. The background from misidentified QCD jets is approximately 
random in sign and so cancels in the t + t~ — r ± r ± combination. To the extent that gg 
and gq production dominates, so does the background where both r's come from xf decay. 
Figure |3| shows an excess of r + r _ pairs ending approximately at the endpoint of the dilepton 
distribution (from ~^ tt) shown in Figure ^[ The structure is less clear due to the energy 
carried off by the neutrinos. 

Even for 5 = there is a substantial rate for real E^rj^ pairs from r pairs produced by the 
decay chain x® ~^ ~ * Xi T+T ~- Two independent chargino decays can also give Et pairs 
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of either sign. Misidentified QCD jets give additional tr^ pairs; the signs of rr pairs from 
such misidentified jets should be random, and this is assumed here. The LFV mixing from 
5^0 produces additional /xr^ pairs, 92% of which have signs that are correctly identified. 
Figure | shows the OS and SS backgrounds and the signal for an assumed 10% branching 
ratio for the direct decay x\ — ¥ Xi^^- This corresponds to 5 = 0.25 (see Figure [[]). Since 
ISAJET assumes lepton flavor conservation, the decay of interest was simulated by finding 
events with two r's consistent with x\ ~^ Xi TT an d with at least one hadronic r decay and 
then replacing the other r with a muon with a probability equal to the assumed branching 
ratio. This is an excellent approximation since, as is indicated above, the mass shift due to 
mixing is very small. Figure || shows the signal and the sign-subtracted background. The 
latter should cancel in the subtracted /jt — er distribution up to statistical fluctuations; the 
LFV signal occurs only in the fir channel. The distribution from the LFV signal has a peak 
at larger values of invariant mass than that of the lepton flavor conserving process as the 
lepton in the latter must arise from the decay of a tau and is therefore softer than that the 
muon from %° — * A* T Xi 

In Figure | for 50 < M tT < 100 GeV, there are 1089 OS and 707 SS events, with equal 
numbers of er and /jt, and there are 518 lepton flavor violating /ir events, 92% of which are 
correctly identified as to sign. Hence in this mass range we expect 

N^rT) = = 0.92(. 5 x 1089 + 518) = 977 
AT( e ± r T ) = .5 x 1089 x 0.92 = 501 

Adding the errors in quadrature, we would then measure 

A^(/i ± r =F ) - V(e ± r T ) = 476 ± 39 

giving a 12.2a excess for 10fb _1 . The statistical 5a limit for 30 fb -1 would be a branching 
ratio of 2.3%, corresponding to 5 ~ 0.1. The signal to background ratio is better in the sign 
subtracted distributions, but the statistical errors are larger. A determination of the branch- 
ing ratio from the observation of an excess of /ir events requires a detailed understanding of 
the systematic uncertainties. The scaling of these results to the design LHC luminosity also 
requires a careful study of the effects of pile up of low-pr hadronic events. It is interesting to 
remark that the sensitivity that we obtain from direct slepton decay is comparable to that 
claimed for a lepton collider operating at 500 GeV center of mass energy for 30 flT 1 [ETI. 



4 Comparison with Rare Decays 

The decay r — ► yq is sensitive to the same lepton flavor violation as the signal considered 
here. The approximate formula of Ref. |^, 

' 5 \ 2 /l00GeV\ 4 



B R(t - 1.1 xlO- 6 ^ ^ ^ 

gives for 5 = 0.1 and M- t = 150 GeV 

BR(t -»• // 7 ) « 1 x 10~ 9 . 
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However this result is sensitive to the details of the mass spectra and mixings; cancellations 
can occur resulting in much smaller rates || pl|. For comparison, the current bound is 
1.1 x 10~ 6 |T^|. The total production of W — > tv at the LHC is about 10 9 events for 
an integrated luminosity of 100 ftT 1 . A study |2(J carried out for the ATLAS detector 
concluded that a 90% confidence level upper limit of 0.6 x 10~ 6 on BR{r — > fij) could be 
reached using W — > tv events for an integrated luminosity of 30 ftT 1 corresponding to three 
years of running at low luminosity. The process is background limited from QED final state 
radiation in the decays W — > \w and W — > tv — > fivv. These contribute approximately 50 
events for 30 fb~ . 



5 Generalization of Results 

We have demonstrated, by using a specific example, how flavor violation in the slepton sector 
can be observed at the LHC. It would be difficult to repeat the analysis for a large number 
of other cases. In order to estimate the generality of the method we have adopted a simpler 
approach. The crucial feature is the production of staus in decays of x\- F° r fixed values 
of A, tan/3 and /i, the parameter space of mo and m\ji was scanned. At each point, all 
masses and branching ratios were computed assuming that there is no slepton flavor mixing. 
If the gluino is lighter than the up and down squarks, then the product branching ratio 
BR(g — > X2 + X) x BR(x2 — ► Tit) is computed, if the gluino is heavier than the up and 
down squarks the combination (BR(u L — > X2 +X) + BR(d L — > \\ + ^0) x BR(x2 —> T i T )/^ 
is calculated. The total supersymmetric cross section is determined at a few points in the 
parameter space and can be approximated by 

a = 1.79 x 10 13 (0.1m + m 1/2 )" 4 ' 8 pb . 

This formula agrees with the cross-section to within 25% for all relavant paramters; the 
dependence on A, tan /3 and \i is small. An approximation of the total production rate for 
f\ from \\ arising from squark and gluino production and decay can then be obtained from 
the product of the cross-section and the combined branching ratio. 

Figures | through |8] show the resulting rates, and Figure [5] shows two sample sets of 
contours of the f\ mass. The light Higgs mass is required to be larger than 103 GeV in 
making these plots; this accounts for the cut off at smaller values of mi/2 for tan/5 = 3. 
Rasing the limit to 107 GeV excludes the region below mi/ 2 = 365 GeV for sgn/i = +. 
The current bound on the Higgs mass for such a small value of tan/3 is approximately 112 



GeV |]I2" |. At this value ra\ii is so large that the production rate is less than 0.1 fb. The 
effect of this bound should be taken cum grano salis as uncomputed higher order corrections 
could modify the relationship between the model parameters and the Higgs mass. The Higgs 
limit has no impact on the plots with tan (3 > 5. The lack of rate at large values of 
is due to the large gluino and squark masses and the consequent small production rate. 
The region where m ^ mi/2 is not accessible as there the decay ~^ TT is n °t allowed 
kinematically. For comparison, the case studied in detail would correspond to a f\ rate of 
~ 5 fb according to Figure [7]. A comparison of this value with the rates shown in Figure |3| 
shows that the combination of acceptance and tau detection efficiency is approximately 
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Figure 6: Contours showing an estimate of the total production rate for f\ as a function of 
m and m 1 / 2 for A = 0, tan/3 = 3, and sgn/i = + (left) and sgn/i = — (right). The contours 
(outer to inner) are at 0.03 0.1 0.3, 1 and 3 pb; the scale is logarithmic. 
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Figure 7: Contours showing an estimate of the total production rate for f\ as a function of 
m and for A = 0, tan (3 = 10, and sgn /i = + (left) and sgn \i = — (right). The contours 
(outer to inner) extend from 0.1 pb to 100 pb and are uniformly spaced on a logarithmic 
scale. 

2%. This is not expected to vary dramatically over the relavent parameter space as the 
kinematics of production and decay are broadly similar. The acceptance will be less near 
the boundary where m^o ~ mf as the resulting taus will be soft. The acceptance is greater 
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Figure 8: Contours showing an estimate of the total production rate for f\ as a function of 
m and mu2 for A = 0, tan/? = 30, and sgn/i = + and sgn = — (right). The contours 
(outer to inner) are at 0.1 0.3, 1, 3 and 10 pb; the scale is logarithmic. 




Figure 9: Contours showing the mass of f\ as a function of m and m\/2 for A = and 
sgn/i = — . Left: tan/? = 3; the contours are 20 GeV apart, from 130 GeV in the bottom 
corner to 610 GeV in the top corner. Right: tan/? = 30; the contours are 25 GeV apart, 
from 75 GeV in the bottom corner to 530 GeV in the top corner. The lower bound in m 1 / 2 
in the left figure results from the requirement > 103 GeV; this bound is less stringent 
for sgn fj, — +. 
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at larger values of ttlq or mi/2 where events are more central and the resulting taus more 
energetic. At a luminosity of 10 33 cm -2 sec -1 a meaningful measurement should be possible 
if the rate is larger than 1 pb. An examination of the plots shows that this covers a sizable 
fraction of parameter space. A detailed quantitative assessment of the sensitive region is not 
straightforward. This is because the background arises from supersymmetry itself and not 
from standard model sources. The background therefore varies across parameter space, so a 
more detailed sensitivity estimate would require repeating the full analysis at each point. 

The identification of 77 decays requires tight cuts on the shower shape and track multi- 
plicity, so it becomes more difficult at high luminosity, 10 34 cm _2 s _1 . For high SUSY masses, 
however, the typically have higher and so are more easily distinguished from QCD jets. 
Extending the analysis described here to high luminosity requires a much more detailed and 
detector-specific study of pileup effects than has been done here. If it proves to be possible, 
then the parameter space where the rate is larger than 0.1 pb should be accessible. This 
covers almost all of the physically interesting region where the decay X2 ~^ T7 ~l * s allowed. 



6 Conclusion 

We have demonstrated that the observation of lepton flavor violation may be possible at 
the LHC We have focussed on the case where the largest violation occurs in the \x — r 
sector. While this is motivated by the atmospheric neutrino problem this violation is harder 
to detect than that occuring in the \x — e sector since the detection of hadronically decaying 
taus is involved. If sleptons are produced in the decays of squarks and gluinos the rates are 
large enough that the resulting sensitivity is greater than can be reached in rare decays of 
the type r — > /ry. This corresponds to the fraction of parameter space where mi/2 > tuq. 
In the SUGRA model the fraction of parameter space where these decays occur with a rate 
large enough for observation is large. 

If sleptons are not produced in these decays, the sensitivity will be much less. A study 
of the possibility of observing direct slepton production for the mass spectrum considered 
(without flavor violation) |TT[. Stringent jet veto cuts are required to reduce the backgrounds 
from other SUSY events. Approximately 50 events from l\ production were found in 30 fb _1 
and a mass sensitivity of 20 GeV obtained. Note that this sensitivity is comparable to the 
mass shifts induced in the l\ sector for 5 = 0.1. A more detailed study is needed before a 
definite conclusion can be made, but it appears that this process might also be more sensitive 
than the r — > p,j decay. 

The importance of 17 decays for SUSY extends beyond the lepton flavor violating signal 
considered here. If tan/3 ^> 1, then the decays xt ~~ * an d xt ~^ can be the only 
allowed two-body decays and so have branching ratios close to unity; then r signatures 
dominate [HJ. For any value of tan/3, 77 signatures can provide information both on the 



splitting between f\ and Ir and on the mixing between tr and 77 [22]. It is therefore 



important to continue a detailed study of 77 signatures for the LHC detectors. 
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